Normal State Correlated Electronic Structure of Iron Pnictides 
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We describe the correlated electronic structure of a prototype Fe-pnictide superconductor, 
SmOi-xFxFeAs, using LDA+DMFT. Strong, multi-orbital electronic correlations generate a low- 
energy pseudogap in the undistorted phase, giving a bad, incoherent metal in qualitative agreement 
with observations. Very good semi-quantitative agreement with the experimental spectral func- 
tions is seen, and interpreted, within a correlated, multi-orbital picture. Our results show that 
Fe-pnictides should be understood as low-carrier density, incoherent metals, in resemblance to the 
underdoped cuprate superconductors. 

PACS numbers: 71.27. +a, 74.25.Jb, 74.70.-b 
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Discovery of high-T c superconductivity (HTSC) in the 
Fe-based pnictides [l[ is the latest among a host of other, 
ill-understood phenomena in d-band oxides. HTSC in 
Fe-pnictides emerges upon doping a bad metal with spin 
density wave (SDW) order at q = (tt, 0). Preliminary 
experiments indicate 0, Q unconventional SC. Existent 
normal state data indicate a "bad metal" without Landau 
Fermi Liquid (FL) quasiparticles at low energy These 
observations in Fe-pnictides are reminiscent of cuprate 
SC. The small carrier density (giving rise to carrier pock- 
ets), along with Uemura scaling from /i-SR [4] similar to 
hole-doped cuprates strongly suggests a SC closer to the 
Bose condensed, rather than a BCS (£ ~ 1000a) limit. 

LDA studies show that SC in Fe-pnictides is associated 
with the Fe-d states hybridized with As-p states: this 
leads to two hole, and two electron- like pockets || . Find- 
ing of van-Hove singularities and peaks at q = (tt, 0) in 
the bare spin susceptibility indicates a SDW state @, 0] , 
found within Hartree-Fock random phase approximation 
(RPA) studies of effective, two- and four-orbital Hub- 
bard models, in apparent agreement with inelastic neu- 
tron scattering results 8] . The observation of quasi- linear 
temperature (T) dependence of the resistivity, a pseudo- 
gap in optics [9(, and a spin gap in NMR [lOj, however, 
is a benchmark feature of the relevance of strong, dy- 
namical spin and charge correlations in the pnictides. 
In cuprates, these are seen in the underdoped state in 
the proximity of a Mott insulator (MI) , suggesting that 
the Fe-pnictides might be closer to a MI than gener- 
ally thought [HI]. Actually, the undoped pnictides show 
an insulator-like resistivity without magnetic order for 
T > 150 K Onset of bad metallic behavior corre- 
lates with a structural (tetragonal-orthorhombic (T-O)) 
distortion at T* ~ 150 K, below which SDW order sets 
in. The small carrier number apparently generated upon 
the structural distortion accords with the observed high 
resistivity, lending further credence to such a view. 

The above suggests that one should study a single Fe- 
As layer with strong electronic correlations to begin with. 
Here, we study the five-orbital Hubbard model within 



the local-density-approximation plus dynamical-mean- 
field-theory (LDA+DMFT) approach, incorporating one- 
electron band structure aspects. Extant LDA+DMFT 
works give either a strongly renormalized FL [l3 | or an 
orbital selective, incoherent, pseudogapped metal 
Apart from the known sensitivity to the value of Jr 
(Hund coupling), the LDA+DMFT spectra show no- 
ticeable qualitative differences at low energy. Are Fe- 
pnictides then strongly renormalized FL metals, or inco- 
herent non-FLs in their normal state? Comparison with 
experimental one-particle spectra should go a long way 
toward resolving this question. We do this in this work. 

Photocmission (PES) and X-ray absorption (XAS) 
studies are a reliable tool to study the correlated elec- 
tronic structure of d- and /-band compounds [3]. 
To date, the only PES experiments have been per- 
formed on SmOi_ x F x FeAs [lj| and LaOi_ x F x FeAs flij] . 
Very recently, XAS has also been performed for 
SmOi_ x F x FeAs pjj ■ Together, they provide additional 
evidence for the "incoherent metal" normal state in Fe- 
pnictides. The PES spectra show a kink at low energy, 
51 = 15 meV [THj], below the T-0 distortion followed 
by SDW order. This kink sharpens with cooling, and 
evolves, apparently smoothly, across T c . Its microscopic 
origin is an enigma. Is it related to the T-0 distor- 
tion, or to the SDW transition? Is it observed only for 
electron-doped systems? Correspondingly, XAS shows 
a well-defined peak at 0.5 eV, and a transfer of weight 
from high energy (2.0 eV) to low energy with F doping, 
a characteristic feature of correlated systems. Answering 
these questions within a correlated electronic structure 
approach provides deeper insight into the underlying cor- 
relations in the non-FL metal phase, aiding in the search 
to identify mechanism(s) of superconductivity itself. 

Starting with the high-T tetragonal structure with lat- 
tice parameters found in Rcf. [18], one-electron band 
structure calculations were performed for SmOFeAs us- 
ing the linear muffin-tin orbitals (LMTO) [T3] scheme 
in the atomic sphere approximation [20j ] . LDA pro- 
vides valuable description of the relevant orbitals on a 
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single-particle microscopic level, but falls short of de- 
scribing dynamical correlations in d- and /-band com- 
pounds. This requires "marrying" LDA to DMFT, 
which enables direct access to the correlated spectral 
functions [2l[. The one-electron part for SmOFeAs is 

H a = Ek.a.ff^W^a.A"," ' where a = x 2 - y 2 ,3z 2 - 
r 2 ,xz,yz,xy label the diagonalized, five d bands. The 
corresponding density-of-states (DOS) (Fig. [1} shows 
that all the five d-bands cross the Fermi energy, Ep , but 
the 3z 2 — r 2 band is almost gapped at Ep. While the xy- 
band has a deep pseudogap at Ep, the x 2 — y 2 ,xz,yz- 
bands have large DOS at Ep. In Fe-pnictides, the d G 
configuration of Fe 2+ dictates that the full, multi-orbital 
(MO) Coulomb interactions must be included. These 
constitute the interaction term, which reads 
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We choose parameters employed by Haule et al. [l2|, 
U = 4.0 eV, U' = U - 2J H = 2.6 eV, and J H = 0.7 eV, 
along with the five LDA bands, and solve H = Ho + Hi nt 
within LDA+DMFT. To solve the MO-DMFT equations, 
we use the MO iterated-perturbation-theory (IPT) as an 
impurity solver. Though not quantitatively exact, it has 
many advantages. It is numerically very efficient, is valid 
at T = 0, in contrast to QMC, and self-energies (E a (u;)) 
can be extracted very easily. These are of particular im- 
portance in a complicated MO situation that occurs in 
Fe-pnictides. 

Given orbital-induced anisotropies in the LDA, strong 
MO correlations renormalize various d-bands in widely 
differing ways. Generically, one expects these to par- 
tially (Mott) localize a subset of d-bands, leading to 
orbitally selective Mott transitions (OSMT), and bad 
metallic states [22[ • This requires strong U,U'. Within 
LDA+DMFT, this orbital selective mechanism involves 
two renormalizations: (a) static (MO Hartree) shifts 
rigidly move various <i-bands relative to each other by 
amounts depending upon their on-site orbital energies 
and occupations, and, more importantly, (6) dynami- 
cal effects of strong U, U', Jh drive large spectral weight 
transfer (SWT) over wide energy scales. Upon small 
changes in bare LDA parameters, large changes in SWT 
lead to OSMT, as well as to incoherent metallic phases 
characteristic of a wide variety of correlated systems. 
With parameters for Fe-pnictides as above, does an 
OSMT occur (HI, or does it not [T3|? What is the origin 
of the observed incoherent metallic behavior? 

To answer these questions, we now turn to our results. 
In sharp contrast to LDA, our LDA+DMFT results show 
drastic modification of the spectral functions. The dy- 
namical correlations lead to dramatic spectral weight re- 
distribution over large energy scales O(5.0) eV. Most 
interestingly, we find no FL quasiparticle signatures in 



FIG. 1: (Color online) Orbital-resolved DOS: LDA and 
LDA+DMFT (U = 4.0 eV, V = 2.6 eV and total band fill- 
ing, n=6.0) for SmOFeAs. The insets show the imaginary part 
of the corresponding self-energies, showing clear evidence of 
non-FL behavior. 



the low-energy spectra; instead, the metallic state is to- 
tally incoherent. In contrast to earlier work (l3j . this 
occurs even without a strict orbital selective Mott local- 
ization, though almost all bands are very close to Mott 
localization. The orbital-resolved self-energies (inset in 
Fig. [T|) clearly reveal this aspect: large damping at Ep 
destroys the FL quasiparticle, giving an incoherent, pseu- 
dogapped, bad metallic state. This incoherent state has 
also been found in previous LDA+DMFT works 12 , l| , 
and, as discussed there, is in qualitative agreement with 
experimental observations. Here, however, we also show 
that our results are in very good semi-quantitative agree- 
ment with key features of the experimental PES and XAS 
spectra, see Fig. [2] 

Analyzing the LDA+DMFT spectra of Fig. [Q we find 
that the 3z 2 — r 2 ,x 2 — y 2 as well as the xz,yz orbitals 
continue to be almost degenerate (the first pair is split 
by 0.06 eV). This is also reflected by the fact that the 
DMFT spectra of 3z 2 — r 2 ,x 2 — y 2 bands show notice- 
able similarities, even as, interestingly, large differences 
between them exist at level of LDA. Of interest are the 
sharp, very low-energy (20 meV below Ep) structures in 



y spectra. In the light of their near- 



the 3z 2 — r",x" 
degeneracy, we ascribe these peaks to the low-energy or- 
bital fluctuations (coupled to charge fluctuations) in this 
two-fold degenerate sector. This is an interesting mani- 
festation of the two- fold 3z 2 — r 2 , x 2 — y 2 orbital degen- 
eracy surviving in Fe-pnictides, and explicitly requires 
strong MO correlations. Using LDA+DMFT, we have 
also estimated the spin state on Fe sites. With parame- 
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FIG. 2: (Color online) Comparison between theoretical 
(LDA+DMFT DOS) and experimental (PES, XAS) spectra 
in the nFL metallic phase of SmOi- x F x FeAs, x = 0.15. No- 
tice the good agreement with experimental PES result up to 
0.5 eV binding energy. PES and XAS data are taken, re- 
spectively, from Refs. [l5l ] and The insets show the In- 
dependence of the experimental angle-integrated PES (left) 
at low binding energies for x = 0.15, and the derivative of 
the total DOS (for n = 6.2). The latter show a sharp peak 
at Q. — 20.0 meV (solid: theory, dashed: experiment), indi- 
cating the kink structure in very good agreement with the 
experimental value of 15 meV below T* . 



ters as above, the quantity < S 2 total > ~ 1.8, slightly 
less than the atomic value of S = 2, and consistent with a 
high-spin state. The observed suppression of the sublat- 
tice magnetization in neutron scattering should therefore 
be attributed to strong geometric frustration in Fe pnic- 
tides, as discussed recently by several authors 0, fill ] - 
This is hitherto an open issue: should Fe-pnictides be 
modelled using S — 1 or S = 2 models [ll|? Our anal- 
ysis resolves this issue in favor of S — 2 modelling. We 
note that this implies strong Jh , and, from d-shell quan- 
tum chemistry, larger U, U' , putting the pnictides into 
the "strongly correlated" class of materials. 

An LDA+DMFT-experiment comparison, carried out 
here for the first time for Fe-pnictides, puts these features 
into deeper perspective. In Fig. ® we compare our results 
to recent PES and XAS results on SmOi_ x F x FeAs with 
x = 0.15. Very good semi-quantitative agreement with 
experiment is visible. In particular, the broad peak at 
-0.3 eV in PES, as well as that at 0.5 eV in XAS, are in 
nice agreement with our calculations. Further, the trans- 
fer of spectral weight from the 2.0 eV feature to that at 
0.5 eV in XAS is also reproduced semiquantitatively in 
the correlated spectra. Repeating the DMFT calcula- 
tions with U = 5.0 eV, as in Ref. 12j, we found worse 



agreement with both PES and XAS data. Most impor- 
tantly, PES reveals a kink at 15 meV upon cooling the 
sample below T* = 150 K, where a T-0 distortion, fol- 
lowed by the q = (it, 0) SDW order, takes place. This 
feature sharpens with decreasing T, and weakens with 
doping, but does not undergo further change across the 
SC T c . This appears to be a more generic feature of the 
electronic structure of Fe-pnictides, as similar evolution 
of the low ener gy p seudogap has also been resolved in 
LaOi- x F x FeAs [3. Analysing our LDA+DMFT spec- 
tra for n = 6.2 (n is the total band filling of the d shell), 
we find a sharp nonanalytic structure in dp total (ui) / duj 
at an energy of 20 meV, implying a kink in the DOS 
at that energy, in very good semi-quantitative agreement 
with PES results. Further, analyzing the orbital-resolved 
DOS, this feature is seen to originate from the 2>z 2 — 
r 2 ,x 2 — y 2 (DMFT) spectra. Since this structure is a con- 



sequence of near two- fold (3z 



-y ) orbital degen- 



eracy, as discussed above, its appearance below T* now 
has an attractive interpretation: it reflects the low en- 
ergy, coupled charge- orbital fluctuations in this pnictide. 
It appears only below T* because the T-0 distortion, in- 
terpreted as a Jahn- Teller instability, occurs at T* [23j, 
and lifts this degeneracy. The itinerant (albeit incoher- 
ent) character of the system suppresses the bare 3z 2 — 
r 2 , x 2 — y 2 splitting to small values (20 meV, as pointed 
out above). Strong orbital fluctuations in this almost 
degenerate orbital sector coupled to one-electron Green 
functions enter the dynamical, second order contributions 
(of the generic form J (uj — uji — cj 2 )g[°' (uj2)G^ (oj — 

LUi)duJiduj2 = J Xab (^i)Ga\uj - uji)dui, with Xab( u ) 
the inter-orbital susceptibility describing dynamical or- 
bital correlations) to the self-energies, and hence show 
up in the LDA+DMFT spectra as sharp, low energy 
peaks in the orbital-resolved spectral functions, as seen 
in our results. We have also computed the LDA+DMFT 
spectra for hole doping (n — 5.8, dot-dashed line in 
Fig. [3]). In contrast to electron doping, no noticeable 
change is observed in the low energy spectra; we predict 
that PES /XAS on hole-doped Fe-pnictides will show this 
in future. 

Since the kink in PES is now related to the T-0 distor- 
tion, rather than the SDW, it should be smeared out for 
T > T* , exactly as observed. Therefore, its survival with- 
out apparent modification across T c is not connected to 
the destruction of SDW order apparently required for SC 
to emerge, //this turns out to be generic for Fe-pnictides, 
it would imply an indirect link, at most, to SC, to the 
extent that it reflects electronic structure changes (viz, 
removal of 3z 2 — r 2 , x 2 — y 2 orbital degeneracy) required 
for the SC instability to emerge from such a normal state. 
"Melting" of the T-0 distortion upon F doping implies 
rapid suppression of this kink feature in our picture: this 
is indeed seen in our LDA+DMFT spectra for n = 6.2 
(orange line) curves. Only a weak remnant of the 15 meV 
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ical correlations are essential to proper understanding 
of the basic physics of Fe-pnictides. As in cuprates, 
the normal state that becomes unstable to unconven- 
tional [ISH SC is not a Fermi liquid. Our LDA+DMFT 
result gives an incoherent, non-FL state with small car- 
rier density, as observed experimentally Very good 
semi-quantitative agreement with extant PES/XAS data 
lend strong credence to our view of Fe-pnictides as 
multi-orbital, strongly correlated materials close to the 
itinerant-localized boundary. 

L.C. thanks A. Ormeci for discussions. L.C. and H.R. 
thank the Emmy-Noether Program of the DFG for sup- 
port. M.S.L. thanks the MPIPKS for financial support. 



FIG. 3: (Color online) Total LDA+DMFT spectral functions 
for SniOi- x F x FeAs. Notice the large dynamical SWT upon 
electron and/or hole doping. LDA result for SmOFeAs is 
shown for comparison. 



kink is resolved at n = 6.2, as shown in the right inset 
of Fig. [21 indicative of surviving uncorrelated distortions 
among an "undoped" fraction of Fe-sites, even as long- 
range distortion melts with doping. 

Based on this detailed theory-experiment agreement, 
we discuss the implications of our work on SC. First, 
since all d bands cross Ep, SC should involve inducing 
the SC gap on all FS sheets. This does not necessar- 
ily conflict with two-band Hubbard model results 241 ] . 
since the multi-band SC proximity effect 25] could op- 



erate here. Once SC pairing occurs in the d xz ^ yz man- 
ifold, such an effect could induce secondary gaps over 
the remaining FS sheets. Second, our results indicate 
that having 2>z 2 — r 2 , x 2 — y 2 orbital degeneracy drives 
a Jahn- Teller T-0 distortion, to the detriment of SC, as 
seen. Finally, our finding of a strongly incoherent "nor- 
mal" state with drastically reduced charge carrier num- 
ber (given the proximity to a Mott insulator in DMFT) is 
consistent with observations [j], 0, 0] indicating a non-FL 
metal with low carrier density SC. Many of these observa- 
tions are reminiscent of those seen in HTSC cuprates [26| 
up to optimal doping, putting the Fe-pnictides into the 
"strongly correlated, unconventional, HTSC" category. 

In conclusion, using a correlated electronic structure 
(LDA+DMFT) approach, we show that strong dynam- 
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